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We apply rotating optical flux lattices to spinor Bose-Einstein condensates. Distinct quantum 
states emerge for fractional ratios of vortex charge density to optical flux density. We exhibit the 
calculated charge-flux states and discuss their topological structure and experimental signatures. 



In the presence of a magnetic field B, an electric charge 
of strength e moving at a velocity v experiences a Lorentz 
force F = ev x B, which causes it to travel in a circu- 
lar trajectory in the plane perpendicular to the magnetic 
field. In strong fields electrons can become trapped by 
the magnetic flux. Quantum mechanically an external 
magnetic field B = VxA couples to the motion of the 
electrons via the magnetic vector potential A correspond- 
ing to the Hamiltonian H = (p — eA) 2 /2m, where p is 
the momentum of an electron with mass to. Integer and 
fractional quantum-Hall effects in a two-dimensional elec- 
tron gas emerge when the number of electric charges per 
magnetic flux quantum, becomes sufficiently small. Typ- 
ically such conditions can be achieved by subjecting the 
sample to a very strong external magnetic field PHI]- 

In scalar superfluids, circulation of the superfluid ve- 
locity is quantized in integer units k = h/m given by 
the Onsager-Feynman rule for quantization of circula- 
tion § v ■ dl = riK, where h is Planck constant and n is 
an integer [5j . This furnishes a topological invariant n in 
the system. When vortices move relative to the flow field 
they are immersed in, they experience a Magnus force, 
which is the hydrodynamic counterpart to the magnetic 
Lorentz force. Within the peculiar analogy between elec- 
tric charges and vortices in hydrodynamics, dating back 
to the works of Helmholtz and Kirchoff, the quantized 
vortices can, in two spatial dimensions, be viewed as the 
analogs of electric charges in a logarithmically interacting 
screened two-dimensional Coulomb gas [6) [7] . 

The connection between superfluid vortices and elec- 
tric charges motivates the search for analogue electro- 
magnetic phenomena in neutral superfluids. Indeed, if 
the fluid is made to rotate at an angular speed fl, the 
Hamiltonian, when expressed in the rotating frame, ac- 
quires an extra term -J] ■ r x p due to the angular 
momentum of the particles. The Hamiltonian H — 
p 2 /2m — ft ■ r x p may be re-expressed in the form 
H = (p„ — kAq) 2 /2m v — mVL 2 r\/2, where m v = 7m 
and p t , = ^/7TO are the effective mass and momentum 
of a vortex defined by the parameter 7, and An = 
— mfi(— ye x + xe y ) is a vector potential. The effect of 
rotation Q on a neutral atom thus appears equivalent to 
the effect experienced by a vortex charge k in a combi- 
nation of a magnetic field Bn = V x An and a scalar 
potential of the centrifugal kind. 

To emulate large effective magnetic fields Bo, neutral 
superfluids must be rotated at high angular speeds ex- 
erting a large centrifugal effect on the superfluid [5HI3]. 
Recently, the creation of synthetic gauge fields utilizing 



the spin-degrees of freedom of neutral atoms [TIMI?] has 
opened up the possibility of generating pure synthetic 
magnetic fields, thereby avoiding the undesired effects 
of the centrifugal scalar potential. Such synthetic gauge 
fields may eventually facilitate the approach to a regime 
of strong atom-atom correlations where the vortex den- 
sity becomes comparable to the atom density, possibly 
allowing the emergence of novel fractional quantum Hall- 
like states of neutral atoms. 

In this paper we investigate the effects of rotating a 
spinor Bose-Einstein condensate using a recently intro- 
duced optical flux lattice [TH [T5] . We study the regime 
where the atoms are weakly interacting and their num- 
ber is much greater than the number of vortices in the 
system. We show that new quantum states emerge in 
this system when the ratio of vortex density to the opti- 
cal flux density acquires rational values. Experimentally, 
such regime is achievable in slowly rotating and optically 
thick atomic clouds. 

We model a gas of ultra-cold atoms cooled below the 
critical temperature for Bose-Einstein condensation in 
the S = 1 hyperfine-spin manifold confined in a parabolic 
potential well T4ra P ( r ) = mujj_(x 2 + y 2 + X z z 2 )/2, where 
luj_ is the transverse harmonic trap frequency and \ z is 
the aspect ratio. In the presence of a field B^, which 
couples to the spin S and rotates at an angular velocity 
fi, we model this system with a Hamiltonian 

H = -h 2 V 2 /2m + y trap (r) - fj, a - S • B«> - ft • L + 

gn(r)+g s J2Ma(r)S„, (1) 

where the total particle density n(r) = ^ \p* a (Y)ip a (r) is 
expressed in terms a three-component spinor order pa- 
rameter = [tpf, ipo, 4>i] T and li a are the corresponding 
chemical potentials and L is the orbital angular momen- 
tum operator. The diagonal operators are expressed as 
scalars. The Cartesian components of the magnetization 
density M a (v) = W{r)S a ^!{?) are obtained as spin-space 
expectation values of the components of the spin-1 oper- 
ator. Depending on the form and origin of B^, the term 
— S • B^ may account for an internal spin-orbit coupling 
or a coupling of a magnetic moment to an external field. 
The repulsive g > particle interactions are modeled us- 
ing the standard contact interaction potential [5] and the 
spin-exchange coupling g s of the atoms is set to be fer- 
romagnetic. The parameters used in our calculations are 
X z = 2, gN/hcu±aQ — 2000, where N is the number of 
atoms in the system, g s — — O.Olg, and do = \Jh/mus^. 
The chemical potentials are m = /ij. ~ \2hio±^ and 
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Mo = l^i + WOhuj±. For the results presented we have 
chosen the scaling clq = 1/im. 

Inspired by the optical flux lattice scheme introduced 
by Cooper and Dalibard jTSJ US], we utilize an optical 
flux lattice of the form 

= g < f > (cos(kx)e x +cos(ky)e v +sin(kx) sin(fcy)e z ), (2) 

where A: is a wavenumber and = 3ftio± determines the 
coupling strength of the optical flux lattice which has 
an areal flux density B^ — l$(k j 'it) 2 . We set the flux 
lattice rotating at an angular speed ft, which controls 
the vortex charge density nucleated in the system. The 
key point here is that it is possible to control Jl and B0 
and hence the number of vortex charges and the num- 
ber of optical flux quanta independently of each other. 
We have shifted the origin of the optical flux lattice by 
7r/2 to set a centre of a flux cell on the axis of rotation. 
That the superposed rotating lattice carries optical flux 
quanta in a form of optical vortices, imprinting geomet- 
ric phases on the condensate, is the crucial distinguishing 
factor between the optical flux lattice scheme [THUS] and 
the conventional, non-topological, optical lattices which 
have previously been applied to nucleate vortex charges 
in Bose-Einstein condensates jTTJ [5D] . Dynamical varia- 
tion of the lattice spacing by a factor of 2.5 has also been 
realized for conventional optical lattices [21]. Similar ex- 
perimental techniques could also be deployed to rotate 
optical flux lattices and to control their flux density. 

We define a charge-flux ratio v = p/q, where p and 
q are the areal density of vortex circulation quanta p = 
£ V N V /A and the optical flux quanta q — i^N^jA, re- 
spectively. In our system the charges are half-quantum 
vortices i v — \ and each flux lattice cell encloses half 
a flux quantum = \ [TBI HH], such that v — N v /N ( j > . 
We extract v from the calculated states by counting the 
number of vortex charges N v inside an area containing 
N<p flux lattice cells. 

Figure 1 (a) exhibits a calculated column density 
n z (x,y) of the condensate integrated along the z-axis 
when an optical flux lattice is applied in a charge neutral 
(nonrotating) system. The corresponding magnetization 
column density |M z (a;,?/)| is shown in Fig. 1 (b), which 
reflects the square lattice structure of the underlying op- 
tical flux lattice. Figure 1 (c) and (d) , respectively, depict 
the condensate and magnetization column densities for a 
synthetically charged (rotating) system in the absence of 
the external optical flux. Circles in Fig. 1 (c) mark the 
locations of the vortex charges whose long-range repul- 
sive vortex-vortex interactions results in their crystalline 
ordering. These vortex charges are half-quantum vortices 
and due to their composite structure, the total atom den- 
sity remains nearly uniform inside the vortex cores. 

When the optical flux lattice is combined with finite 
vortex charge density, charge-flux correlations develop 
in the system resulting in quantum phases which fea- 
ture self-assembled topological structures. Charge-flux 
commensurability is achieved when the ratio of the vor- 
tex charge density n v , which may be estimated by the 
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FIG. 1. (Color online) Flux-only and charge-only states. Con- 
densate column density (a) and magnetization density (b) for 
a charge neutral system in the presence of an optical flux. 
Condensate column density (c) and magnetization density (d) 
for a synthetically charged system in the absence of optical 
flux. The circles in (c) mark the locations of the half-quantum 
vortex charges. The rotating condensate (c) has a larger di- 
ameter compared to that of the nonrotating cloud (a) due to 
the centrifugal effect. 



Feynman's rule n v — 2VL/k to the optical flux density 
B(j, is a rational number 1/i — n v /B^. For our sys- 
tem the corresponding commensurable wavenumbers are 
fcj = 2iTyJi£l/ n. On varying k and/or Q, commensuration 
energy related to inter-vortex interaction and topological 
pinning of vortex charges by the optical flux quanta is 
anticipated to drive the system into new ground states. 
Effects of commensurability between a vortex lattice and 
a non-topological optical lattice have been studied previ- 
ously in both weakly interacting and strongly correlated 
systems and the results are reviewed in |12j . 

Figure 2 shows numerically calculated states with dif- 
ferent values of v obtained for fixed ft = 0.7a; ± and by 
varying the optical flux density B^. We choose to fix 
the value of ft since the flux density is proportional to 
the square of k whereas the vortex density is linear in 
ft. Thereby we also avoid having to deal with the large 
variation in condensate and vortex densities, which are 
strongly dependent on ft. The first column of Fig. 2 
shows an integrated column density n z (x,y) of the con- 
densate calculated for different optical flux densities. The 
second and third columns, respectively, display the den- 
sity of condensate atoms in two of the hyperfine sub- 
states, which have practically equal populations. The 
third hyperfine state has been lifted in energy with re- 
spect to the other two states and is not shown due to 
its negligible population. Experimentally, this could be 
achieved using a quadratic Zeeman effect. The third col- 
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FIG. 2. (Color online) Synthetic charge-flux states with 
different fractional values of v. Condensate column den- 
sity (first column) and corresponding spin-up (second col- 
umn) and spin-down (third column) component densities and 
magnetization column density (third column) and its Fourier 
transformation (last column) corresponding to different ratios 
v marked on each row (a)-(d). Circles in second and third 
column frames mark the locations of the half-quantum vortex 
charges. Two Frenkel-type defects are highlighted by a cross 
in the fourth frame of row (c). Optical flux densities in (a) - 
(d) are B4, = {0.10,0.20,0.46 and 0.62} (jim)" 2 , respectively. 
The optical flux lattice rotation frequency is Q — 0.7uj± for 
all of these states. 



umn shows the magnetization column density \M z (x,y)\ 
and the last row displays its Fourier transform. In Fig. 
2 (a) each cell of the flux lattice traps two half-quantum 
vortices realizing a v = 2 state while in the v = 1 state 
(b) every flux lattice site is occupied by a single half- 
quantum vortex charge. Generically states with integer 
v are composed of interlaced lattices of i^-fold quantized 
vortices. In the half-filling state (c) every second row of 
the flux lattice is free of vortex charges. The topology of 
the fractional v — 1/3 state is shown in the last row (d). 

By inspecting the states shown in Fig. 2 we extract a 
trial spinor ^P(r) = /(r)[^> Q (r),^g(r)] T , where T denotes 
a transpose and the function /(r) is a smooth envelope 
determining the finite system size. Here the spinor com- 
ponents are Vv(r) = Yl v ( z — z v) p " Ylj e - '^' ^ 2X \ where 
j and v respectively enumerate the flux cells and the cells 
occupied by charges in the state a. For v < l,p v = 1 and 
p v = v for integer values of v . The length A is a vari- 
ational parameter and the index v labels the complex 
coordinates z v = x v + iy v of the vortex charges. 

These quantum liquids support macroscopic topologi- 
cal fractional-vortex excitations which emerge as impu- 
rities such as interstitials and vacancies. A pair of such 
Frenkel-type defects are highlighted in Fig. 2 (c) where 
the two half-quantum vortices shown appear rotated 90 




FIG. 3. (Color online) Ratio \ jv (a) and an orbital angular 
momentum per particle L z /N as functions of the optical flux 
density £?</,. The calculated data points are joined by solid 
lines. The vertical dashed lines are plotted at commensurate 
flux densities 13$ with ki = 2-K^JiQ./ hi where i = 1/2, 1, 2, 
and 3. The horizontal base line in (b) is drawn at the value 
corresponding to the angular momentum in the absence of the 
optical flux lattice shown in Fig. 1 



degrees about their centre-of-mass — half way through ex- 
changing their equilibrium places. Experimentally, con- 
trolled vortex braiding comprising of adiabatically mov- 
ing the vortices [22] around each other would require 
addressability of the optical flux lattice at a single site 
level. In addition to the macroscopic topological defects, 
these quantum liquids possess a low energy spectrum of 
microscopic Bogoliubov-de Gennes collective quasiparti- 
cle modes. Of particular interest is the level structure 
of those quasiparticle modes associated with the optical 
flux quanta, vortex motion and the edge localized high 
angular momentum surfon states. 

Figure 3 (a) shows a measurement of the ratio of 
the optical flux density to the vortex density \jv as a 
function of the optical flux density B$. Figure. 3 (b) 
shows a measurement of the expectation value of the z- 
componcnt of the orbital angular momentum per particle 
L z /N as a function of the optical flux density B$. The 
vertical dashed lines in both subfigures indicate the com- 
mensurate flux densities estimated using the Feynman's 
rule for the average vortex density. The commensurate 
wavenumbers fc, seem to be slightly overestimated by this 
method indicating that the number of vortices in this sys- 
tem is not large enough for it to accurately mimic solid 
body rotation. The Feynman's rule has also been ob- 
served to overestimate the vortex number in experiments 
[llj . For commensurable ratios of 1/v we observe min- 
ima in Fig. 3(b) coincident with the plateaux in Fig. 3(a). 
The horizontal line in Fig. 3(b) corresponds to the angu- 
lar momentum of the flux- free state shown in Fig. 1. 

The obtained results can be understood in terms of 
interaction of the vortex charges with the optical flux 
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quanta. At commensuration, the ratio of the average 
vortex charge density to the average optical flux density 
is a rational number, corresponding to the plateaux in 
Fig. 3(a)). The optical flux lattice is then rotating at 
an optimal angular speed with respect to its flux density 
such that the condensate superflow and the transport 
of vortices are in balance. Increasing the optical flux 
density introduces compressive stress to the vortex liquid 
forcing it to increase its density, which is tantamount to 
creating excitations and radial vortex currents. Higher 
vortex density manifests itself as an increase in the orbital 
angular momentum. 

Since the angular speed of the vortices is restricted by 
the fixed rotation frequency of the flux lattice, the super- 
flow can not be in balance with the flow of the vortices for 
incommensurate charge-flux ratios. The relative motion 
between vortices and condensate flow gives rise to mutual 
friction [2 6) in the presence of excitations. The vortices 
resist compression by creating vortex currents which re- 
store the balance between vortex density and flux density. 
This facilitates transitions between quantum states with 
different lattice symmetries shown in Fig. 2. 

The staircase behavior seen in Fig. 3(a) and the cor- 
responding angular momentum fluctuations shown in 
Fig. 3(b) result from the alternating compression and re- 
laxation of the vortex charge density. The compressibil- 
ity of the background gas of condensate atoms due to the 
gapless sound wave spectrum does not exclude the pos- 
sibility of coexistence of an incompressible vortex liquid 
since the part of the excitation spectrum which couples 
to the motion of vortices (Kelvin- Tkachenko modes) may 
be gapped and decoupled from the compressional sound 
waves. 

The predicted charge-flux states shown in Fig. 2. can 



be prepared experimentally by applying an optical flux 
lattice [HI [19] to spinor Bose-Einstein condensates and 
rotating it at a fixed angular speed. Such rotating drive 
has previously been achieved using non-topological opti- 
cal lattices [TTJ [2U] . Transitions between different charge- 
flux states can be observed by varying the wavenumber 
of the optical flux lattice [3T]. Subsequently the conden- 
sate density can be imaged using standard absorption or 
phase-contrast imaging techniques, either in-situ or af- 
ter a time-of-flight and the different spin-components can 
be resolved individually by separating them using Stcrn- 
Gerlach fields. The ratio v of vortex charge density to 
optical flux density can be extracted from the obtained 
images. The magnetization structure which clearly re- 
veals the distinct lattice symmetries of the charge-flux 
states could be directly measured using magnetization- 
sensitive phase contrast imaging [23l [24]. The orbital 
angular momentum can be inferred from spectroscopic 
measurements of the collective excitation spectrum of 
the condensate [25] ■ Resonance spectroscopy can also 
be used to explore the low-lying excitation spectrum of 
these quantum liquids. 

In conclusion, we have computationally applied a ro- 
tating optical flux lattice to spinor Bose-Einstein conden- 
sates and have found distinct quantum states and order- 
ing of the nucleated fractional-charge vortices at rational 
ratios of the number of vortex charges to the number of 
optical flux quantum. These charge-flux states could be 
observed in slowly rotating and optically thick cold atom 
clouds using current experimental technology. It will be 
particularly interesting to probe the excitation spectra to 
reveal the gapped states and zero modes of these quan- 
tum liquids. 
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